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Abstract 

This article presents a direct measurement of the growth angle during the growth 
of a cylindrical 2" silicon crystal using a radio- frequency heated floating zone 
process. From the high-resolution pictures taken during the process, this growth 
angle was evaluated to be 11° ± 2°. Furthermore, the free surface of the melt 
was modeled using the Laplace- Young equation. This model has to include the 
electromagnetic pressure calculated by the surface ring currents approximation. 
The results were compared to the experimental free surface derived from video 
frames. It could be shown that the calculated free surface will only fit the 
experimentally determined one if the right growth angle is considered. 
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1. Introduction 

The floating zone (FZ) technique is used to grow silicon single crystals from 
polycrystalline feed rods (see ref.[l] for details). FZ silicon is mainly used for 
high-power electronics and is standing out against crucible-grown Cz silicon 
through the low oxygen content and high electron lifetime. A detailed descrip- 
tion is found in ref.[2]. 

Investigating the free surface of the melt, one advantage of the FZ process is 
the good visibility. One can take photographs of nearly the full contour of the 
melt zone. For optically heated FZ furnaces, some investigations about the free 
surface can be found in ref.[3] for silicon and GaAs or in ref.[4] for sapphire. In 
our case, videos are normally taken to control the growth process by measuring 
zone height and diameter of the crystal. In this paper, we used video frames to 
compare the numerical calculation of the free surface, based on the paper of Lie 
et al. [5], with the experiments. 

The second task of the paper is the verification of the growth angle. There are 
some evaluations of the angle, e.g. in references [6, 3, 7]. For the FZ process, 
this was described for processes in mirror furnaces. In our investigations we 
use radio-frequency (RF) heated FZ furnaces. Observing the FZ growth of 



Email address: wuenscher@ikz-berlin.de (Michael Wiinscher) 
link to the original document: 10.1016/j Jcrysgro.2010.11.101 



Preprint submitted to Elsevier 



January 12, 2013 



silicon in such a furnace, the presence of a growth angle could be doubted. 
Differences to observations in mirror furnaces could have been explained by the 
electromagnetic pressure. However, careful analysis of the captured free surface 
as well as the high-resolution photographs, enabled the determination of the 
growth angle. 



2. Model for the free surface of the melt 



The FZ process is nearly axisymmetric and, therefore, the pressure difference 
along the free surface of the melt was treated that way. The Laplace- Young 
equation 1 is solved on the geometry shown in fig. 1. 

AP-pgz + F n = 7 K , (1) 

with the pressure difference to surrounding gas AP, the radial/vertical coordi- 
nates r/z, the electromagnetic pressure F n , the mean curvature the surface 
tension 7, the density of the melt p and the gravitational acceleration g. The 
equation is rewritten as a system of differential equations in order to be solved 
with numerical methods [5]. 



feed rod 




1 - solidification interface 

2 - free surface 

3 - crystalisation interface 

4 - symmetry axis 

5 - open melt front 
ITP - internal triple point 
ETP - external triple point 



Figure 1: Scheme of the FZ process, melt zone with the dimensionless coordinate system 
(grey) and the second point (R S ec, Z se c), if the solution fits. 
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The dimensionless coordinates R/Z are scaled with the crystal radius Rc and 
the system has two material-dependent parameters, the Bond number Bz and 
the electrical Bond number Bei depending on the electrical angular frequency 
lj, the permeability \i and the applied voltage to the induction coil Vb- 

This system of equations is solved by the Runge-Kutta method with a dis- 
cretization step of 0.002 for the dimensionless arc length. We solve the system 
as an initial value problem with a start angle /?o at position (0., 1.) until cur- 
vature / angle at the end position are reached and the calculation is stopped. 
Typically the end angle is chosen to be zero, this means, that the curve ends 
vertically. 

The equation has one free parameter Pq, which is fixed by choosing a second 
point where the solution has to go through. This is accomplished by an opti- 
mization algorithm. From an initial solution, where Pq is chosen to be not too 
far from the final solution, the distance to the present second point (R se c, Z sec ) 
is calculated. This point is determined by selecting the first point on the solu- 
tion, which fits the radial coordinate R sec . If the solution does not have such a 
point, the last point of the solution is selected. Therefore the last term of P^ ew 
directs the solution towards the goal point. This term is zero for the last steps 
of convergence. 

r>New _ tdOU % — Z sec R — R sec 

The growth angle is, as an initial condition for the solution, an important pa- 
rameter for calculating the free surface. With a smaller start angle the solution 
has a lower melt height, which can be seen in figure 2. 




Figure 2: Free surface calculation for different start angles (0° and 11°) without electromag- 
netic pressure. Same neck diameter of 10mm with vertical ending for a 2" FZ silicon process. 



3. Model for Electromagnetic Pressure Calculation 

The FZ process is heated by a RF induction coil at about 3MHz. From this 
coil a current is induced into the melt. As a side effect, there is the Lorentz 
force pressing against the free surface, which is taken into account with the third 
term of equation 4 on the right-hand side. In this paragraph, we will show a 
way to calculate the current. 
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With respect to the skin depth 5 = y (5cu,3MHz ~ 0.04mm) the cur- 

rent can be treated as surface current. For axisymmetry, it is further assumed 
that this surface current is modeled by ring currents [8] around the z-axis in 
cylindrical coordinates. With these assumptions, the dimensionless surface cur- 
rent Jq is calculated as follows: 

A e - C ° S ^ [ Je{s)\j ^[K{rn) - E(m)]ds (6) 




m = . (7) 

where R/Z are the coordinates of the ring current (source), which produces a 
magnetic field quantified by the vector potential at the coordinates r/z. All 
vector potentials are superimposed by integration along the surfaces of the con- 
ducting materials and give the equation for the system vector potential Aq with 
the integration along the arc length s. In cylindrical coordinates this potential 
contains the functions K(m) and E{m), the complete elliptic integrals of first 
and second kind. Along the surfaces of crystal, melt and feed rod, where no 
voltage is applied, the vector potential Aq is zero. Only along the induction 
coil, the high-frequency voltage gives a vector potential Aq = CQ ^.^ . 

Together with the boundary conditions, we achieved an integral equation 
for the unknown surface current Jo(s), which can be numerically solved by the 
boundary elements method. Therefore, the geometry is split into linear line 
segments numbered by n, where the dimensionless current Jg n is assumed to 
be constant and equation 6, only depending on the geometry, is solved with a 
10-point Gaussian quadrature at the center of the line elements (r^, Zi), yielding 
the coefficients A-„ . 



^segments 




, induction coil 



/ A in J Gn = < \ \ (9) 

7^1 1° > else 

The given system of linear equations, numbered by z, is solved with a standard 
direct numerical solver for the unknown currents Jg n . A detailed description 
can be found in the article of Lie et. al. [5]. Our program is written in python 
using the scipy [9] and numpy libraries as well as wx Widgets for the interface. 
A speed up of about 5-10 times is reached by using Fortran for calculating the 
integral. 

By solving this linear system, we get the discrete form of the surface current 
Je(s). One numerical difficulty is to map this discrete solution to match the val- 
ues required by the numerical solver of the Laplace- Young equation. Therefore 
we transform the numerical values from the cylindrical coordinate system into 
a parametrization along the arc length s. These values are then interpolated 
by a spline function of third order, which is used as input for the Runge-Kutta 
solver of the Laplace- Young equation in the 2. paragraph. We toggle between 
both solvers until we do not have to change Po for a given error of 0.001 in 
dimensionless units. 
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4. Comparison of numerical results with experiments 

To achieve pictures of the experiments, we used our camera system for mea- 
suring the diameter of the crystal during the growth of a 2" crystal. These 
black and white picture have a resolution of 720 x 576 pixels what is 5.83 pixel 
per mm. For a horizontally aligned camera, the inner end of the free surface is 
hidden by the induction coil. Therefore we cannot use the end of the surface as 
second point (R seci Z sec ). Instead, we used an edge detection algorithm, which 
finds the point at 95% of the visible height of the melt surface. It is not the 
last visible point because there are reflections of the melt at the induction coil, 
which make it impossible to find an accurate edge. The solution depends on the 
material parameter density g = 2530 -H and the surface tension 7 = 0.783 — 
[10]. 



(a) p = 0° without EMP (b) p = 11° without EMP (c) p = 11° with EMP 

Figure 3: Comparison of the numerical solution with the captured picture of the free surface. 
Electromagnetic pressure (EMP) (Vb = 500V) needs to be included to get a good fit with the 
experiment. 

The first fits were made without including the electromagnetic pressure. Two 
solutions for initial angles of (a) 0° and (b) 11° are shown in figure 3. Obviously, 
these solutions do not fit well. The deviation is especially large in the middle 
of the curve and even bigger for the solution with the right start angle. If the 
solution includes electromagnetic pressure of the right strength of 500V, the 
solution fits very well for the given resolution of the pictures (fig. 3(c)). 

As shown in the paper of Lie et al., the highest strength of the electromag- 
netic pressure is generated in the middle of the curve. This is also shown in our 
calculations. So the highest change of the free surface is achieved in this region. 
The electromagnetic pressure pushes the solution against the free surface. Of 
course, then one has to find a new Po, else the second point would not lie on 
the curve anymore. 

With the same electromagnetic pressure, the start angle of 11° shows the 
best fit to the experiment (see figure 4). So, the growth angle known from the 
literature could be verified. For the lowest angle (fig. 4 a) the calculated solution 
is below the real melt surface and for the highest angle (fig. 4 c) it is over it. 
It is very interesting to mention that the growth angle is not directly visible 
in the pictures because of the bad resolution, but this angle has a significant 
impact on the whole shape of the free surface. Therefore the growth angle is an 
important parameter to find the coincident free surface. 

Nevertheless, in the applied way, this method has no better accuracy com- 
pared to a direct measurement of the growth angle as described later. In addition 
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(a) po = 0° 



(b) po = 11° 



(c) po = 20° 



Figure 4: Different start angles with the same voltage Vq = 500V for the electromagnetic 
pressure. Best equality with the experiment is given for a start angle of 11°, which is equal 
to the growth angle, (published in [11]) 

to Po, the voltage for calculating the electromagnetic pressures is a second pa- 
rameter, which has to be fitted. The dependency is shown in figure 5. Therefore 
one gets a parameter range for angle and voltage, where both parameters are not 
independent. Other uncertainties are geometry and position of the induction 
coil, which have to be modeled with a sufficient accuracy. 




(a) Po = 9° V = 400V (b) p = 11° V = 500V (c) p = 13° V = 600V 



Figure 5: Best solution with different start angles po and strengths Vo of the electromagnetic 
pressure. Solution b) gives the best result. 



5. Measurement of the growth angle from high-resolution photographs 

As already mentioned, the growth angle is not directly visible in the cylin- 
drical growth of the FZ crystals. Additionally, the small ripple of the crystal 
surface, which is caused by rotation, sometimes enhances this impression. Even 
looking at low-resolution video camera pictures does not show a visible angle, 
but evaluating high-resolution pictures from the experiment, this impression 
changes. 

We have used a Canon EOS 350 D with a Pentax 135mm lens mounted 
with three distance rings of together 50mm and an adapter Hama M42 to the 
camera. The equipment was put onto a tripod in front of the window in the 
growth chamber, which is around 350mm away from the growing crystal. It 
was separately tested that this equipment has a resolution of about 66 pixel 
per mm. Taking photos of the FZ process, it is not possible to use automatic 
focusing because there are many reflections. For taking good photographs we 
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are focusing with a lower aperture, which has a lower depth of focus than a 
higher aperture. The maximum aperture was then used to take the photo at an 
exposure time of l/160s. To reduce the impact of the lens distortion, we placed 
the area of interest in the middle of the photo. 

As in figure 6(a) the crystal radius is just reducing, so the outer shape is not 
vertical anymore. Therefore, the growth angle must not be measured against 
the vertical, but against the tangents of the crystal surface starting from the 
three-phase line. The edges get more visible from the Sobel gradient of the 
picture and two tangents are fitted (see 6(b)). One tangent goes from the phase 
transition point along the free surface and the other one along the crystal. The 
angle between both tangents is the growth angle. 




(a) Picture 



>1 

(b) Gradient 



Figure 6: (a) Picture of the growing 2" silicon crystal (b) Gradient of the zoomed picture with 
tangents to measure the growth angle (11.1°) 



Because of the dynamics of the growth process, this angle is not constant. 
There are vibrations caused by rotation as well as jerks by the movement of 
crystal and feed rod. Also the three-phase line has some grooves at the facets, 
which break the rotational symmetry and will cause a small error in evaluation 
if the photo is taken nearby. During the experiment, we took several pictures 
similar to fig. 6(a) with a delay of some seconds, so we had 4-5 pictures per 
rotation. This was done at the beginning of the cylindrical growth and repeated 
after ten minutes. The growth angles measured from these photos are given in 
the graphs of fig. 7. The mean growth angle is 11° ± 2° for the first series 
and 11.5° ± 1.9° for the last one. Theses values are in good agreement with 
the literature. In one work about the solidification of a wafer heated with an 
electron beam, the growth angle was found to be 11° [6]. A newer work [7] 
of the solidification of a melt on top of a crystal found an equal result of 12°. 
For mirror furnaces, as mentioned before, Tegetmeier [3] got a slightly different 
value of 8.6° for silicon. 

For a single evaluation of the growth angle, two angles need to be extracted 
from the gradient pictures. The reproducibility for that is better than 0.5°, 
therefore, the error for a single growth angle measurement is better than 1.0°. 
This contribution to the mean value is small and so we use the standard devia- 
tion as a good error estimate for our experiment. There are also some systematic 
errors which were tried to reduce. See experimental description above for details. 
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(a) Begin of cylindrical growth 
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(b) End of cylindrical growth 



Figure 7: Measured growth angles with mean value (horizontal line) at two different stages of 
the growing crystal 



6. Conclusion 

We have achieved a good agreement of the simple model using the Laplace- 
Young equation together with the ring current approximation and the captured 
free surface of the experiment. For our fit we had to include the growth angle 
as initial condition as well as the effect of the electromagnetic pressure for the 
best results. Based on our good results we are looking forward to use our model 
to improve our control system for automatic crystal growth. It is also planned 
to increase the resolution of the camera system to get more measurement points 
along the free surface for a more precise comparison. 

We have shown that the growth angle can be detected in RF-heated floating 
zone crystal growth and that it should be measured against the changing crystal 
surface. Although this angle is hardly seen in our experiments, this angle was 
evaluated to be 11° ±2°, which is in the range of the literature value of 11°. More 
accurate measurements could be achieved by a higher resolution of the camera, 
maybe with an automated system for the analyses, by reducing the movement 
of the melt by optimizing the growth parameters like pull speed, rotation rate 
or crystal radius, and by changing the lens to get a bigger depth of focus. 
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